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ABSTRACT 


A  numerical  comparison  of  the  exact  solution, 
the  uniform  asymptotic  solution,  and  the  non- 
uniform  asymptotic  solution  of  a  field  that  produces 
a  circular  caustic  is  performed.  It  is  seen  that  the 
uniform  asymptotic  solution  is  accurate  even  at  the 
caustic  for  moderate  values  of  ka.  Curves  and 
data  are  given  as  a  function  of  ka  and  distances 
from  the  caustic.  The  asymptotic  expressions  are 
derived  both  from  the  exact  solution  and  from  the 
physics  of  the  circular  caustic. 
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SECTION  I 


DERIVATION 

4 


The  uniform  asymptotic  expansion  developed  in  Lewis'  paper*  is  valid 
in  the  neighborhood  of  a  caustic,  and  it  reduces  to  the  usual  geometrical 
optics  representation  away  from  the  caustic.  In  Appendix  I,  these  ex¬ 
pansions  are  calculated  for  the  circular  caustic.  Herein,  an  exact  so¬ 
lution  of 

(f  +  k2)  n  =  0  ,  (1) 

whose  asymptotic  expansions  agree  with  that  of  Appendix  I,  enables  us  to 
obtain  a  numerical  comparison  of  an  exact  solution  and  its  uniform  and  non- 
uniform  approximations. 

An  exact  solution  of  (1)  in  two  dimensions,  where  polar  coordinates 
are  used,  is  the  function 

ikab  T  .  /0, 

At  =  e  Jka(kr)-  (2) 

Here  J  (z)  is  the  Bessel  function  of  order  n.  We  obtain  a  uniform  asymp- 
n 

totic  expansion  of  (2)  by  using  the  following  formula,  given  in  Asymptotic 

t 

Approximations  by  H.  Jeffreys: 


*  *  R.  M.  Lewis,  Uniform  Representation  of  Geometrical  Optics  Fields 

Near  a  Caustic,  The  MITRE  Corporation,  ESD-TR-65-404,  Bedford, 
Mass.  ,  May  1966. 

*  t 

H.  Jeffreys,  Asymptotic  Approximations,  England,  Oxford  University 
Press,  1962,  p.  84,  Equation  (36). 
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In  (3),  Ai  is  the  Airy  function  defined  and  discussed  in  Jeffreys’  work*  and 
in  Appendix  II.  The  formula  is  valid  for  large  n. 

In  order  to  apply  (3)  to  (2),  one  sets 


z  -  sech  i u  =  (cos  v)  ,  v  =  lju 


(4) 


Thus, 


3/2 


2 1 3  £  =  -iv  -  tanh  (-iv)  =  i  (tan  v-v)=  i^Jz2  -1  -  cos  *  1  / z^  , 


2  2  2 
tanh  pt  =  1-sech  =  1-z  , 


(5) 


and 


J  (nzW 
n 


4  £ 


-.1/4 


1-Z 


-1/3  /  2/3 

n  '  Ai 


("2/30- 


(6) 


Here 


-£  =  3/2 


>jz2  -1  -  COS  1  (l/z^ 


2/3 


(7) 


in  agreement  with  (5).  Let  n  =  ka,  z  =  r/a  and  apply  (6)  to  (2).  The  result  is 


ikad 

4£ 

1/4  / 

V  .2/3  " 

M~e 

1 

>-$  1 

< 

P 

to 

1 _ 

(ka)  /  Ai 

* 

^ _ 1 

(8) 


or 


ik0  ..  (  .2/3  \ 
e  g  Ai  (^-k  /  p  J 


(9) 


*See  page  1,  second  reference. 
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Here 


0 


+ 


and 


0  =  a0  ,  (10) 

p  =  -a2/3=  (3/2  xt  )2/3  ,  2/3  p3^2  =  ^  =  \!r2  -a2  -  acos  1  (a/r)  ,  (11) 

1/4 


g  =  (ka) 


-1/3 


-4£ 


21  2  n 
r  /a  -1 


„  2/3  -.1/4 

4k  P  n 


[k2(r2-a2)J 


(12) 


The  non-uniform  expansion  is  obtained  from  (9)  by  using  the  asymptotic 
expansion 

Ai  (-£  )  ~  £  1^4sin  (2/3  £  3^“  +  7r/4^;  -  7r/3  <  arg  £  <  tt/3;  |e  |  — (13) 


of  the  Airy  function.*  Let 

£  =  k 2/3p  .  Then  2/3  £  3^2  =  k4<  and 


1  ikab 

^  I 6 

^  7T 


-,1/4 


2  2  2 
k  (r  -a  ) 


sin 


+  7r/4  j  , 


(14) 


or 


2  7Ti 


(r2-‘ 


1-1/4 


)eik(ap+t)+i  tM  _eik(a^)-i  Wi)  (15) 


This  result  may  be  interpreted  in  terms  of  geometrical  optics  with 
the  aid  of  Figure  1 . 

Two  rays,  each  tangent  to  the  circular  caustic  of  radius  a,  intersect 
at  the  point  X  whose  polar  coordinates  are  r,  0.  The  two  phases  associ¬ 
ated  with  these  rays  are 


*  See  page  1,  second  reference,  p.  33,  Equation  (30). 
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Figure 


Circular  Caustic  Corres¬ 
ponding  to  Equation  (15) 


(16) 


where 


=  a( <t>  +  a)  =  a<?  +  a  cos 


-1 


(17) 


We  see  from  (11)  that 


ap  +  4*  =  aP  +  a  cos 


-1 


T  =  S. 


(18) 


hence,  (15)  can  be  written  in  the  form 


i7r/4 

e 

— — 

M  2  7rk 


iks-  ,  iks  ,  -  i7r/2  , 
z  e  +  z+  e  + 


(19) 


% 


4 


4 


where 


z  =  z 


(20) 


+ 


In  (20)  ,  the  amplitude  and  the  phase  of  the  incoming  and  outgoing  waves  are 
clearly  exhibited.  The  amplitudes  become  infinite  at  the  caustic  (t=0),  and 


-iv/2 


the  phase-shift  factor  e 


multiplies  the  outgoing  wave.  The  wave 


fronts  are  the  spirals  S+  =  constant. 

Here,  the  uniform  expansion  (9-12)  has  been  obtained  directly  from 
the  exact  solution  (2).  In  general,  of  course,  the  exact  solution  is  not  known. 
Then,  it  is  necessary  to  use  the  method  developed  in  Lewis'  report*  to  ob¬ 
tain  the  uniform  expansion.  In  Appendix  I,  that  method  is  applied  to  this 
problem,  and  a  uniform  expansion  is  obtained  that  agrees  exactly  with  (9-12). 

Numerical  computations  and  graphs  of  the  exact,  uniform,  and  non- 

uniform  results  will  not  be  obtained.  For  this  purpose,  0=0  may  be  taken, 

since  a  common  factor  of  all  three  results.  It  is  noted  that  ka, 

2/3 

kr,  k'i'  ,  and  k  p  are  dimensionless,  therefore  let  k=l.  The  exact, 
uniform,  and  non-uniform  results  are  denoted  by  p  ,  p  and  p  .  Then, 

1  Z  tj 


(21) 


(22) 


and 


1/4 


(23) 


# 


See  page  1,  first  reference. 
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Here 


and 


(24) 

(25) 


<* 
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SECTION  n 


NUMERICAL  RESULTS 

These  three  functions  were  programmed  for  several  values  of  ka  as 
a  function  of  r/a.  Figure  2  shows  these  curves  for  ka  =  2,  5,  and  10.  For 
these  values  of  ka,  curves  of  the  exact  and  uniform  solutions  are  indis¬ 
tinguishable,  and  only  curves  of  #x  and  /x  are  shown.  To  indicate  the 

J-  O 

errors  in  the  uniform  approximations,  a  deviation  curve  is  plotted  for 
ka  =  1  and  ka  =  2  in  Figure  3.  Figure  4  is  the  same  as  Figure  2  except 
that  ka  =  l/2,  l/4,  l/8,  and  l/l6. 

For  these  values  of  ka,  u  has  diverged  sufficiently  from  /jt  so 

^  1 

that  three  curves  can  be  distinguished.  From  these  curves  the  following 
conclusions  are  clear.  The  non-uniform  approximation  is  very  accurate 


* 


i 


# 


Figure  2.  Plot  of  Field  versus  Normalized  Distance  for  Various  ka 
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♦ 


V 


Figure  3.  Plot  of  Error  versus  Relative  Distance  for  Two  Values  of  ka 


i 


Figure  4.  Plot  of  Field  versus  Normalized  Distance  for  Various  ka 
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« 


* 


8 


for  sufficiently  large  r/a  and  all  ka.  The  uniform  approximation  is  good 
for  almost  all  r/a  even  for  quite  small  ka.  To  indicate  the  actual  errors 
in  the  uniform  approximation,  the  maximum  full-scale  percentage  error  for 
various  ka  is  indicated  in  Table  I.  This  is  defined  as 

maximum  |  n  -y.  1 100 

u  X 

maximum 

From  the  curves  of  Figure  3,  it  can  be  seen  that  the  deviation  from  the  ex¬ 
act  solution  at  a  given  point  may  be  considerably  smaller  than  this  value. 

Examination  of  the  field  on  the  caustic  using  (21),  (22),  and  (23)  leads 
to  the  exact  solution  =  J^a(ka),  while  the  uniform  approximation  is  inde¬ 
terminate  and  the  non-uniform  approximation  is  singular.  In  Appendix  III, 
an  expression  is  developed  for  the  uniform  approximation  (26). 


Table  I 

Percentage  Error  for  Various  ka 


* 


* 


ka 

Percentage  Error 

1/64 

87.5 

1/32 

57.  6 

1/16 

31.  6 

1/8 

17.  3 

1/4 

8.  05 

1/2 

3.3 

1 

1.  7 

2 

0.45 

3 

0.  27 

5 

0.  18 

10 

0.  082 

20 

0.  037 

50 

0.  017 

9 


(26) 


21/3 

"a1/3 


'-2l/3£' 

a1'3 


In  Figure  5,  the  actual  percentage  error  at  the  caustic  as  a  function  of  ka* 
is  plotted.  Since  the  largest  value  of  the  exact  solution  is  not  at  the  caustic, 
these  errors  are  slightly  larger  than  those  of  Table  I.  As  ka  becomes 
large,  the  error  goes  to  zero.  As  ka  goes  to  zero,  p  becomes  singular. 

Li 

At  small  ka  the  exact  solution  approaches  unity,  and  hence  the  relative 
error,  (m2  - >  approaches  the  difference. 


* 


4 


Figure  5.  Plot  of  Error  on  the  Caustic  as  a  Function  of  ka 

_  « 

*  This  curve  was  suggested  by  D.  Ludwig,  Courant  Institute,  New  York,  N.  Y. 
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APPENDIX  I 


THE  CIRCULAR  CAUSTIC 


j 


Herein,  aspects  of  the  application  of  various  representations  of  geo¬ 
metric  optics  fields  near  a  caustic  shall  be  considered.  The  standard 
asymptotic  expansion  and  a  representation  due  to  R.  M  Lewis*  shall  also 
be  considered.  To  be  specific,  the  two-dimensional  field  that  produces  a 
circular  caustic  was  chosen. 

The  wave  front  that  produces  a  circular  caustic  is  a  spiral  determined 
by  a  circle  and  constant  length  string.  In  Figure  6,  we  have  the  circular 
caustic  of  radius  a.  The  points  1*  and  2’  lie  on  the  spiral  wave  front. 

They  are  determined  by  the  condition 


(27) 


where  L  is  a  constant  and  p is  the  distance  on  the  wave  front  to  the 
caustic;  i.  e. ,  p^  is  the  magnitude  of  the  line  segment  11*  . 


The  field  at  any  point  off  the  caustic  is  determined  by  the  two  rays 
that  pass  through  the  point.  In  Figure  2,  the  two  indicated  rays  pass 
through  the  point  P.  The  distance  from  the  point  i  on  the  wave  front  to 


P  is  designated  by  r  .  The  distances  from  P  to  1  and  2  are  equal  and 
are 


(28) 


where  rka  is  the  radial  coordinate  of  the  point  P.  We  note  that 


(29) 


(30) 


* 


See  page  1,  first  reference. 
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Figure  6.  Circular  Caustic 
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f 


* 


J 


> 


or  that  one  ray  through  p  has  touched  the  caustic.  The  angular  location  of 
points  1  and  2  are  found  from 

0  =  0  +  cos  1  a/r  ,  (31) 

where  0  is  the  angular  coordinates  of  P.  Using  (30),  (27),  and  (28)  on 
(29),  the  following  is  obtained: 

=  L-a(0  -  cos  1  a/r)-  ^Jr"2  -  a2  ;  (32) 

=  L-a  (0  -  cos  ^  a/r)+  'jr2  -  a2  .  (33) 


The  standard  asymptotic  formula*  for  the  field  at  p  is 


u  = 


elk(*2+T2)>  (34) 


where  the  subscript  refers  to  the  ray.  This  expression  differs  from  that  in 
Lewis'  report*  since  this  is  a  two-dimensional  problem  where  p  is  used 
instead  of  the  curvature  k  .  The  term  associated  with  the  second  ray 
involves  taking  the  square  root  of  a  negative  quantity.  This  indicates  the 
phase  shift  of  71/2  at  the  caustic.  This  formula  is  not  valid  at  the  caustic 
since  it  is  singular  at  such  points.  Using  the  previous  results  for  the 
circular  caustic,  relation  (34)  becomes 


^ik^s^  +  L  -  a  0  +  a  cos 


+  L 


-  a0  -  a  cos  1  a/r  +  ^r 


* 

See  page  1,  first  reference,  Equation  (93). 


(35) 
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For  computation,  the  following  assumptions  are  convenient, 


♦ 


Z.  e 


iksi 


=  Z„  e 


iksr 


=  k 


1/2 


-  a  cos  1  a/rj 


\  l/2 


(36) 


0  =  0.  (37) 

Equation  (37)  requires  that  the  field  is  constant  along  the  wave  front,  and 
the  results  must  be  normalized  so  that  this  constant  is  unity.  That  0  =  0 
follows  from  a  reasonable  choice  of  coordinate  system. 

To  simplify  the  result,  L  can  now  be  chosen  such  that 


* 


ikL  -i7r/4  , 

e  e  =  1  . 


(38) 


Then  (35)  reduces  to 

elkL  e+1  tt/4  .  -ln,  ,  vr~2  2* ,  ,.\ 

u  =  2i  ^2^  2 - 2^1 /4  sin  v^a  cos  /r  "  *  Mr  "  a  +  7r/4 J  * 

and  Equation  (39)  is  rewritten 


(39) 


u  = 


k1'2  (r2  -  J) 


■2^1/4  sin  (ka  cos  *  a/r  -  k\r^  -  a^  +  71-/4J. 


(40) 


Lewis*  gives  a  uniform  asymptotic  expansion 


ikb-i 
Z,  e  1 

N 

CD 

5? 

& 

_ 1 

Ai  |-£2  | 

1 

2 

1  -  rip. 

1/2  + 

~  :  ri/2 

1  -  T  /p.  ' 

f  1 

1/1 

2/2 

i  Ai’ 


ikb-, 


Z^  e 


r,  ikb 

V  2 


-v4/2 


\ 


(41) 


*  See  page  1,  first  reference,  Equations  (68)  and  (69). 
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y 


where 


-J 


¥ 


£  =  (3  k /4)^3(T2  -  T^1^3 


(42) 


and 


h  T!  +  T2 

b.  =  s.  + - 

1  l  2 


(43) 


Using  Equations  (32),  (33),  (36-38),  and  (41), 


^  _  2TT1!2  (3W2)  (\  r2  -  a2  -  a  cos  1  a/r  ) 


1/6 


.  l/2  /  2  2\  1/4 

k  \r  -  a  ^ 


Note  the  disappearance  of  the  derivative  of  the  Airy  function  because 
of  the  symmetry  of  the  caustics.  If  the  asymptotic  expansion  of  the  Airy 
function*  is  used  in  (44),  Equation  (40)  is  recovered. 

The  comparison  of  Equations  (40)  and  (44)  with  Equations  (22)  and  (23) 
shows  agreement  except  for  a  numerical  constant  of  'J  l/27r  .  Thus,  a  check 
is  provided  on  the  validity  of  the  manipulations.  This  factor  arises  from 
the  arbitrary  choice  of  amplitude  for  the  wave. 


w 


* 

See  page  1,  second  reference,  p.  33,  Equation  (30). 
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APPENDIX  II 


THE  AIRY  FUNCTION 


The  Airy  function  is  defined  (e.  g.  ,  in  Jeffreys'  work)  as 

e  t,"1/^t3  dt  , 

L. 


Al(z)  ‘  sr  j 


(45) 


where  is  the  contour  in  the  t-plane  shown  in  Figure  7. 

The  transformation  t  =  is  is  introduced.  Then  for  real  z  the  trans¬ 
formed  contour  can  be  shifted  to  the  real  axis,  and 


Figure  7.  Contour  Relevant  to  the 
Airy  Function 


*  See  page  1,  second  reference,  Equation  (1). 
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00 


Ai(z) 


—  oo 


z  real 


(46) 


is  obtained. 

An  alternative  representation  is 

oo 

Ai(z)  =  ~  \  cos(sz+s^3)  ds  ;  z  real 


(47) 


v 
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APPENDIX  HI 


THE  UNIFORM  ASYMPTOTIC  FIELD  AT  THE  CAUSTIC* 


A  quick  examination  of  Equation  (44)  or,  equivalently,  (22)  does  not 
indicate  that  the  uniform  asymptotic  expansion  is  finite  at  the  caustic. 

In  Lewis'  report^,  it  was  shown  that  this  was  generally  true.  Here,  this 
shall  be  shown  directly  from  Equation  (22).  The  resultant  formula  will  also 
give  a  simple  expression  near  the  caustic.  Let 


r  =  a  (1  +  e)  , 

where  e  is  a  small  positive  quantity.  Then 


2  2  2 
Nr  -a  =  aN2e+e  ~a\2e  , 


(48) 


(49) 


and  from  relation  (24) 


t  \  2  2 

*  -  Nr  -  a  _  -1  ,  ,  v 

-  -  cos  (a/r) 


(50) 


Taking  the  cosine  of  both  sides  and  expanding  the  power  series  to  second 
order  in  e  and  first  order  in  4*  , 


= 


_  2 {2  g3^2 


3  1/2 

a 


(51) 


is  obtained.  For  the  field  near  the  caustic,  the  following  can  now  be 
obtained. 


[-©'I 


(52) 


*We  wish  to  thank  J.  D.  R.  Kramer  of  The  MITRE  Corporation,  Bedford, 
Mass,  for  this  derivation. 

t  See  page  1,  first  reference,  Equation  (6). 
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